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Abstract 

Exactly solvable models provide an opportunity to study different as- 
pects of reduced quantum dynamics in detail. We consider the reduced 
dynamics of a single spin in finite XX and XY spin 1/2 chains. First 
we introduce a general expression describing the evolution of the reduced 
density matrix. This expression proves to be tractable when the com- 
bined closed system (i.e. open system plus environment) is integrable. 
Then we focus on comparing decoherence and thermalization timescales 
in the XX chain. We find that for a single spin these timescales are 
comparable, in contrast to what should be expected for a macroscopic 
body. This indicates that the process of quantum relaxation of a system 
with few accessible states can not be separated in two distinct stages - 
decoherence and thermalization. Finally, we turn to finite-size effects in 
the time evolution of a single spin in the XY chain. We observe three 
consecutive stages of the evolution: regular evolution, partial revivals, ir- 
regular (apparently chaotic) evolution. The duration of the regular stage 
is proportional to the number of spins in the chain. We observe a "quiet 
and cold period" in the end of the regular stage, which breaks up abruptly 
at some threshold time. 



1 Introduction 

In the present contribution several results and observations are presented which 
concern the dynamics of exactly solvable models and especially one-dimensional 
spin chains. The integrability allows to obtain rigorous results, however we 
believe that the applicability of general conclusions based on these results is not 
restricted by integrable models only. The paper is organized as follows. In the 
next section the quantum Liouville equation is rewritten in the form which may 
be tractable when the model is exactly solvable. In section 3 a sketch of the XY 
model is given. In section 4 a comparison of decoherence and thermalization 
timescales in the isotropic XY (i.e. XX) model is carried out. Section 5 is 
devoted to the dynamical regular-to-chaotic transition in the XY chain. In 
section 6 we focus on the quiet and cold period in evolution of a spin in the 
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XY chain which takes place in the end of the regular evolution. The details of 
calculations are presented in two appendices. 

2 Reduced dynamics of a spin: general expres- 
sion 

Consider a quantum system (described by a Hilbert space S), which interacts 
with a quantum environment or bath (described by a Hilbert space B). The 
composite system with Hilbert space "H = S <E> B is considered to be closed. Its 
density matrix evolves according to 

p(t) = cxp(-iHt)p{0) cxp(iHt). (1) 

The Hamiltonian may be decomposed as 

H = H s + H B + H SB , (2) 

where H s and H B are self-Hamiltonians of the system and the bath corre- 
spondingly, and H SB is an interaction Hamiltonian. Here and in what follows 
the usage of superscripts and subscripts S, B is believed to be self-explanatory. 
The spectral decomposition of the total Hamiltonian reads 

d 

H=J2En\*n)(*n\, (3) 

71=1 

where ty n are the eigenvectors and d = dim'H = dgdg. 
The time evolution of p s {t) reads 

p s w = EE (* n |p(o)i* m ) e - i(£ "- £ ™ )t pL, (4) 

n— 1 m— 1 

where the dg x dg matrices 

Pnm = trB|*n)(*m| (5) 

are introduced (not to be confused with matrix elements!). Evidently, these 
matrices encode the dynamics of the open system S, while the matrix elements 
( v E , „|p(0)|5 , m ) describe the initial conditions. 

In the case when the system S is represented by a single spin, its reduced 
density matrix p s (t) and matrices pf m may be parameterized by polarization 
vectors: 

p s {t) = (1 + p(i)<r)/2, p s nm = {S nm + Pnm <r)/2. (6) 

The polarization vector of a spin p(t) belongs to a unit sphere which is known 
as the Bloch sphere. The length of the polarization vector equals 1 for a pure 
state and is less than 1 for a mixed state. Its time evolution is given by 

d d 

PW -EE (*n|p(0)|* m >e-^- £ '" )t Pn m . (7) 

n—1 m—1 
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As a rule this expression is not very useful in practical calculations (however, it 
is used to prove a necessary condition for thermalization of a spin coupled to a 
bath, see ref. pQ). However, when the model is integrable, one has a chance to 
calculate p nm explicitly. If this is done and if the initial condition is of some 
simple form (e.g., the bath is in a thermal state), then eq.Q may appear to be 
tractable and allow to study the reduced evolution of a spin in detail. We will 
demonstrate this is the case of integrable one-dimensional spin chain - the XY 
chain. Before we turn to this specific model let us give an example of initial 
conditions that simplify the problem greatly. These conditions correspond to 
the uncorrelated system and environment, when the system is in arbitrary pure 
state, and the environment is in thermal state with infinite temperature: 

p(0) = 2- N (l + Pq <t)®1 b . (8) 

In this case eq.© reads 

d d 

p(f) =2- N J2Yl (P°> P™)e- <(B "-^ )t p„ m . (9) 

n— 1 m— 1 



3 XY model on a circle 

Consider a chain of N coupled spins 1/2 with the following Hamiltonian: 

N , N 

H=l $> + tK**,. x + (1 - TK^x) + \ £ oj. (10) 

3=1 J'=l 

Here the index N + 1 is identified with 1, and N is supposed to be even. The 
following parameters enter the Hamiltonian: 

• ft — coupling constant, 

• 7 - anisotropy parameter, 

• h - magnetic field. 

We consider the case of weak and moderate coupling: h > k. 

For the purposes of the present study we regard the first spin as the system 
S, and other (N — 1) spins - as the environment B. 

Finite XY model is " almost diagonalizable" through the sequential Jordan- 
Wigner, Fourier and Bogolyubov transformations (see e.g. [2] [5]). An important 
property of the XY Hamiltonian is that it commutes with the parity operator 

[Il,H]=0. (11) 
The diagonalization leads to the "almost free fermion form" of the Hamiltonian: 

H = p° dd E q { c + q c q -\) + p cv E E M c «-\)> ( 12 ) 

q£X odd q£X ov 

where 

N N N N 1 N 3 N 1 

X dd - { — -+1,-^+2,...,—}, X ev = {-—+-,-—+-,...,— --}, (13) 
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Figure 1: Patterns of the spin polarization evolution p(i). Left: what may 
be expected for the relaxation of a spin in the analogy to the relaxation of 
macroscopic bodies ("artist's view"). Right: real relaxation of a spin in the XX 
chain for h — 10k (plot). The points correspond to the values of p(i) from t = 
to t = 4k" 1 with the step 0.01k" 1 . The initial value reads p(0) = (0.6,0,0.8). 
See the text for other details. 

{c q , q € X 0( jd} and {c q , q <E X cv } are two sets of fermion operators (note, 
however, that the operators from different sets do not satisfy anticommutation 
relations, see ea. (|4"5|) '). P odd and P ev are parity projectors, 

p cv = (i + n)/2, p odd = (i-n)/2 (w) 

and fermion energy is defined as 

E q = ^{h - KCOs(2vrg/iV)) 2 + ( 7 k sin(27rq/^)) 2 . (15) 

The details of the XY spin chain diagonalization are presented in the ap- 
pendix A. 

4 Decoherence versus thermalization timescales 

Decoherence is a process of the dynamical vanishing of non-diagonal entries of 
the reduced density matrix p (t) written in the pointer basis. Thermalization 
is a process of the energy exchange between the system and the environment 
leading to thermal equilibrium. What is the relation between the timescales of 
these processes? Evidently, for macroscopic bodies the decoherence timescale 
is much less than the thermalization timescale. Namely, decoherence timescale 
should be microscopic, in order to ensure the fast quantum-to-classical transition 
as it is required by our everyday experience. In contrast, the thermalization 
timescale is usually macroscopic. In fact, this hierarchy of the timescales is the 
basis for the very distinction between the decoherence and the thermalization: 
the former emerges as a first short stage of quantum relaxation, and the latter 
- as a succeeding long stage. 
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Is the above reasoning applicable to small quantum systems, in particular, to 
a single spin 1/2? Very often one implicitly assumes that it is. However, below 
we present a counterexample: the decoherence and thermalization timescales of 
a spin in the XY spin chain appear to be of the same order. 

The easiest way to show this is to consider a special case with zero anisotropy 
(7 = 0, XX chain) and small coupling (k <C h). The latter condition ensures 
that the pointer basis is the eigen basis of the self-Hamiltonian \a\ [I] in the 
sense that p s (t) is approximately diagonal in this basis for almost all times [5] . 
Let us take the initial state to be 



which corresponds to the environment at zero temperature and arbitrary po- 
larization po of the first spin. We are interested in the evolution of the central 
spin polarization p(i). In the analogy to the relaxation in macroworld one could 
expect the following relaxation pattern (see the left picture on Fig. 1). First 
p x and p y would vanish abruptly (or, equivalently, p s {t) would become approx- 
imately diagonal in the pointer basis {| I ti)})> which would constitute 
decoherence. In contrast, p z would decrease gradually revealing the energy flow 
from the first spin (i.e. system) to cold environment. This pattern would allow 
to unambiguously separate decoherence from thermalization, as is evident from 
Fig. 1, left. 

However the described above pattern is not realized in reality. Instead de- 
coherence and thermalization proceed simultaneously, see Fig. 1, right. To 
be more exact, the process of relaxation can not be separated in time in two 
different stages - decoherence and thermalization. This can be seen from the 
approximate analytical solution which is valid for < t < N / n [B] : 



where Jo is the Bessel function. 

One can see that the relaxation of all components of p(t) occurs at a com- 
mon timescale ~ In a system which consists of a single spin there is no 
large dimensionless parameter which could be responsible for the hierarchy of 
timescales. This reasoning suggests that the absence of such hierarchy is a gen- 
eral property of quantum systems with small number of accessible states rather 
than a peculiar feature attributed to integrability of the XY chain. However 
further work is necessary to show this with more rigor. 

5 Finite-size effects in XY chain 

Integrable spin chains are often studied in the N —¥ 00 approximation. Al- 
though this allows to use the powerful machinery of Quantum Field Theory, 
some interesting features of the evolution may disappear in this approximation. 
The present section is devoted to one of such features - the regular-to-chaotic 
transition during the evolution of p(t). A similar effect in a somewhat different 
setting is observed (but not discussed in detail) in [7]. 

In what follows we take the initial state © which corresponds to the environ- 
ment at T = 00. As was already discussed such choice of the initial conditions 



p(0) = (1 + Po<ti)/2 ® I UU ... Iw)<| 2 | 3 ... 4-jv I, 



(16) 



p x = Jo{nt)(pQ cos ht — Pq sin ht), 
p y = Jo(nt)(pQ sinht +Pq cosht), 
P z = -l + (l+pg)J?(/st), 



(17) 



5 




Figure 2: Evolution of p z (t) for 7 = 0, h = 5, re = 1, N = 100. 

simplifies the problem. However one should be cautious when interpreting the 
infinite temperature state 2~ N+1 1 B . Indeed, consider a thermal state with some 
high temperature T. It is well approximated by 2~ N+1 1 B whenever T ^> Nh but 
not just T ^> h. Evidently the former condition (in contrast to the latter one) 
becomes more stringent with the growth of N and appears to be impracticable 
in the limit N — > 00. However this subtlety is not important in our discussion: 
the qualitative picture of the regular-to-chaotic transition does not depend on 
the temperature of the environment. 

We focus our attention on the evolution of the z-component of polarization. 
Using the explicit expressions for the eigenstates of the XY model one gets from 
eq.© the following exact formula (the details of calculations see in appendix B): 

P z (t) = \po(A 2 odd + A 2 CV + B 2 odd + B 2 CV ), (18) 

where 

Aev(t) = N- 1 J2 cos E q t, A odd {t) = N- 1 J2 cos E q t, 

B ev (t) = N- 1 cos 0q sin Eqt, B odd (t) = N^ 1 £ cos 6 q sin E q t, 
qex ov qex odd 

(19) 

h — ncos{2iiq IN) , , 

cos 9 q = v q/ ' . (20) 

Eg 

The plot of p z {t) for 7 = (XX chain), h = 5, re = 1, N = 100 is presented 
on Fig. 2. One can distinguish three successive stages of evolution: 

1. < nt < N : regular evolution. At this stage p z (t) is well described by 
the approximate expression [5] 

p*(i)~^J 2 (rei). (21) 
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Figure 3: Evolution of p z (t) for various anisotropics and magnetic fields, n — 
1, N= 100. 

2. N < Kt < few iV : partial revivals. 

3. nt > few iV : apparently chaotic evolution. 

As the duration of the regular stage is proportional to N, it is the only stage 
which may be catched by the N — > oo approximation. Evidently this stage 
ends up when the perturbation initially localized in the first site of the chain 
propagates through the circle and returns back to the first site. However the 
abruptness of the deviation of p z (t) from the regular expression (12"T1) at some 
threshold time calls for additional explanation, which will be the purpose of the 
future work. 

The second stage has a well-defined starting point, but its termination point 
is not well-defined. One can clearly see one partial revival at N < nt < 2N in 
Fig. 2, but the succeeding evolution can not be characterized that unambigu- 
ously. In fact the second stage gradually changes to the third one - apparently 
chaotic. Note that due to the Poincare recurrence theorem the revivals should 
interrupt the chaotic evolution at long times. 

On Fig. 3 one can see the plots for p z (t) for various anisotropies and mag- 
netic fields. Evidently large anisotropy suppresses the oscillations, while large 
magnetic field amplifies them. 

6 Quiet and cold period 

The asymptotic behavior of the Bessel function is given by 



Jo (a;) ~ t/^/ttx cos(x — 7r/4). 



(22) 
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Therefore, according to eq. (f2"Tj) , during the regular stage the smoothed behavior 
of p z (t) in the XX chain is described by 



™ = St (23) 
At the same time the long-time average of p z (t) is given by 

p z ~ 2p z /N. (24) 

Note that at times N(2ttk)^ 1 < t < Nn" 1 the smoothed p z (i) is Zess than the 
average value p z . This means that during this time interval the energy of the 
central spin is less than its long-time average energy. Also the oscillations of 
p z {t) are suppressed at this time. Thus we observe the quiet and cold period 
which precedes the onset of chaotic evolution. Such " calm before the storm" may 
be of practical importance in the context of quantum information processing. 



7 Acknowledgements 

I am grateful to E. Safonov for useful discussions. I acknowledge the partial 
support from grants NSh-4172. 2010.2, RFBR-11-02-00778, RFBR-10-02-01398 
and from the Ministry of Education and Science of the Russian Federation under 
contracts N^NP 02.740.11.5158, 02.740.11.0239. 



Appendix A: diagonalization of finite XY spin 
chain 

The diagonalization of the XY spin chain was studied in many works starting 
from the pioneering paper [5] . The early papers focused on the thermodynamic 
N — > oo limit. The finite size XY chain is discussed in later works, see e.g. [5] [3] . 
Here we present the details of diagonalization for the sake of completeness. 

A.l H in terms of aj 

We define the operators op in a usual way, 

These operators are neither Bose nor Fermi operators: 

<^7 + <W = 1 > (*) 

= °"n CT / for 3 ^ n - ( 27 ) 

The following simple equalities prove to be useful: 

a z = 2a + a- - 1 = -2(j-<j+ + 1 (28) 
a z a+ = -a+<r z = a+ , o z o~ = -a-<j z = -a~ (29) 
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The Hamiltonian may be rewritten in terms of aj 1 as follows: 



H = H a + H 7 + Hh 



with 



N 

j=l 

N 

K l V~V + + - - \ 

H ~t = > ■^■+1 + ^^+1). 

AT 

H h = h^af(ij -Nh/2. 

A. 2 Jordan- Wigner transformation 

Dchne the operators 

n 

Evidently, II jy coincides with the parity operator II defined in section 3. 
Define Fermi operators aj, a+ as follows 

a 7 = a j~ U i-i = U i-^ a J^ a t = a j" n i- 1 = U 3- 1(J t 

This implies 

aj = ajlij-i, 0+ = afllj-i, 
{«!.<} = <W, {a„,a+} = {a m ,a-} = 0, 



*5=2a+aj-l = -2aja+ + l. 

The Hamiltonian takes the form (note that now the ordering of a^, 
important; also note the change of the total sign): 

N 

H ° = ~ 2 E( a i~°7+i + a i+i°7) ~ + n )( a w«r + 



AT 

ff,, = /i^ataj -iWi/2. 
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A. 3 Fourier transformation 

Define for arbitrary real q 



„;tt/4 n 



i-E^Je-^ ^' N a- nl (42) 
J+e^^^MI/V, (43) 

V n— 1 

Then 

{6+ 6+} = {6fc,6"} = 0, (44) 

2 1 _ g27ri(fc-g) 
i h k > 6 9 } = TV 1 - e 2wi(k-q)/N ■ ( 45 ) 

In particular, if one takes 

AT AT AT 

9 = -2+1. -y+ 2 '-'T ^d) (46) 

N 1 ^ 3 AT 1 

9 = y + 2' -t + 2'-'T _ 2 (Xov) (47) 

then the set of b~ is the set of Fermi annihilation operators. 
The Hamiltonian may be written in terms of 6+ as follows: 

H = H odd P odd + H CV P CV (48) 



with 
and 



p° dd = (l-n)/2, P cv = (l + n)/2, (49) 



#o° dd(CV) = -« E 008(2^)6+6", (50) 

(cv) 

ffodd(ev) = l K7 £ sin(27r g /AT) (6+6t g + 6l g 6-), (51) 

<?6 -X' odd (cv) 

H odd(c V ) =h b+b--Nh/2. (52) 

<?e-fodd (cv) 

A. 4 Bogolyubov transformation 

Dchnc the following quantities 



T q = ^Ksm(2nq/N), e q = h - KCOs(2nq/N), E q = + T\, (53) 

9 q = -arcsm(T g /E q ). (54) 
The Bogolyubov transformation reads 

9 

Cq = cos b~ - sin 6t r (55) 
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The odd and even parts of the Hamiltonian take the form 

#odd(ev) = £ E q {c+c--\). (56) 

This completes the diagonalization. 
A. 5 Eigenstates 

Let us first construct the (pseudo)vacuum states with respect to the annihilation 
operators c~, i.e. the states |vac) dd, |vac) cv which satisfy 

c~|vac) odd (cv ) =0 V 9 € X odd (cv) . (57) 

This states may be written as 

|vac) odd = ^oddcI Ar/2+1 c: Ar/2+2 ...c w/2 | tt ••• t}> (58) 

|vac) ov = ^cvC_ N /2+l/2 C -N/2+?,/2--- C N/2-l/2\ ^ 1") ' (^9) 

where K dd and H ov are normalization constants. This constants arc finite, which 
may be seen from the equalities 

KU-4|c: Ar/2+1 c: JV/2+2 ... c - /2 |tt-t)l = J] cos y^°> ( 6 °) 

q£X odd 

KU-4|c: Ar/2+1/2 c: JV/2+3/2 ...c- /2 _ 1/2 |tt-t)l= n cos ^"^°- ( 61 ) 

qex ov 

Note that |vac) cv is indeed an eigenstate of the Hamiltonian, while |vac) dd is 
not. 

All the eigenstates of the Hamiltonian are obtained from the vacuum states 
by applying the creation operators c+ . To create the odd number of fermions one 
should use q € X Q dd and |vac) dd, while to create the even number of fermions 
one should use q G X cv and |vac) ov . 

Evidently one can enumerate all the eigenstates of the Hamiltonian by the 
multiindcxes 

Qm = {qi, q2, 9m}, 0<M<iV, (62) 

with the ordering q\ < q 2 < ... < qu- Then the eigenstate with M fermions 
read 

\Q M ) = 6+6+...6+ M |vac) odd ( ev) (63) 

with gi,g2, ■••,9m € X oAd ( ov ) when M is odd (even). 

For our purposes we need only the matrix elements between the states with 
the same parity, therefore we use the notation |vac) without subscripts in what 
follows. 
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Appendix B: calculation of p z (t) 

Here we outline the basic steps in obtaining the main result of section 5 - the 
formula for p z (t) in case of the infinite initial temperature of the environment, 
i.e. eqs.([T5]).(P]l. 

The polarizations which enter the general expression ([7]) now should be writ- 
ten as 

P a Q M Q a =&MK\QM). (64) 

In particular, 

■PqmQm = (QmI V«r - MQm), (65) 
p QmQm = ^M\ a t + h\Qm), P V Qm q £i = ^QaMi - (4\Qm)- (66) 
Evidently, p z - ^ only when M — M is even, while p' y ~ ^ only when 

QmQm QmQm 

M — M is odd. Taking this into account one obtains from eq.© 

Qm,Qm 

where E(Q M ) = £ E q- 

The polarizations p z ~ are nonzero in three cases: 



1. 

PQmQm 



v(Qm,P)cos6 p , (68) 



N 
pex 



2. 



where X M = X odA ( cv ) if M is odd (even), and t](Qm,p) = 1 if P S Qm 
and —1 otherwise. 



} M = K M -i U {g}, Q M = K M -i U {q}, q,q<£ K M -i, q ^ q ■ 
Qa/ = Qm-2 U {k} u k,k<£Q M -2, k^k: 

\Pn n \ = \PA n | = -| 1 Sin ^ ~ ^ ] (70) 



Using the above formulae one obtains eas. (|T8|) . (|19[) from ca. (|67[) after straight- 
forward (although somewhat bulky) calculations. 
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